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We have performed calculations of the size of the frequency shift induced by a static electric field 
on the clock transition frequencies of the hyperfine splitting in Yb + , Rb, Cs, Ba + , and Hg + . The 
calculations are used to find the frequency shifts due to blackbody radiation which are needed for 
accurate frequency measurements and improvements of the limits on variation of the fine structure 
constant, a. Our result for Cs (8u/E 2 = -2.26(2) x 10~ 10 Hz/(V/m) 2 ) is in good agreement with 
early measurements and ab initio calculations. We present arguments against recent claims that the 
actual value might be smaller. The difference (~ 10%) is due to the contribution of the continuum 
spectrum in the sum over intermediate states. 

PACS numbers: 32.60.+i,31.30.Gs,31.25.Eb 



I. INTRODUCTION 

Atomic clocks are now important for both practical ap- 
plications and fundamental physics. The hyperfine struc- 
ture (hfs) transition of the ground state of 133 Cs serves 
as a primary frequency standard providing the definition 
of a metric second. Many similar hfs transitions in other 
atoms and ions are used or under consideration for use 
as secondary microwave frequency standards. Most fre- 
quency standards (atomic clocks) operate at room tem- 
perature. However, the exact definition of a metric sec- 
ond corresponds to the frequency of the transition mea- 
sured at zero temperature. This means that readings 
from atomic clocks should be corrected to account for the 
effect of black body radiation (see e.g. The value of 
this effect can be found from measurements or calcula- 
tions. There are many experimental H IE El IE IE an d 
theoretical El IE IE EE El El El works studying the 
effects of black body radiation on microwave frequency 
standards. However, the situation is far from being satis- 
factory. There is disagreement among the different works 
for cesium which will be discussed in more detail below. 
On the other hand, the data for other atoms and ions is 
very poor or absent. 

In the present work we perform accurate calculations 
of the shift due to blackbody radiation in the hyper- 
fine transitions in Yb + , Rb, Cs, Ba + , and Hg + . These 
transitions play key roles in microwave frequency stan- 
dards. They are also of interest to physicists consider- 
ing experiments to measure a variation in the labora- 
tory (see e.g. 0, 0, 0, [l7]). An experiment is cur- 
rently planned utilizing Yb + |18|. No accurate calcula- 
tions or measurements of the radiation shift for this ion 
have been performed previously. However, calculations 
or measurements are available for other atoms and ions 
El II H II IE IE S HIE EE El El El We compare our 
results with these values. 

Due to the importance of cesium as a primary fre- 
quency standard we have performed a more detailed 
study of this atom. This was the subject of a separate 
short paper EE]- In the present paper we give more de- 
tails of the calculations while also discussing other atoms 



and ions. 

There is some disagreement on the actual value of 
the radiation frequency shift for cesium. Early measure- 
ments d H IE] an d ab initio calculations [E EE] support 
a value which is close to —2.2 x 10~ 10 Hz/(V/m) 2 while 
more recent measurements [E and semiempirical cal- 
culations |E fill IT^ | claim that the actual number might 
be about 10% smaller. While we cannot comment on 
the experimental results, the source of disagreement be- 
tween theoretical values seems to be in the continuum 
spectrum. Ab initio calculations (including the present 
work) include the contribution of the continuum spec- 
trum into the summation over the complete set of states 
while semiempirical calculations do not. We demonstrate 
that adding the contribution of the continuum spectrum 
to where it was missed brings all theoretical results into 
good agreement with each other and with early measure- 
ments. 



II. THEORY 

Blackbody radiation creates a temperature dependent 
electric field, described by the Planck radiation law 



7r exp(a;/fcsT) — 1 



(1) 



This leads to the following expression for the average 
electric field radiated by a black body at temperature T: 



(E 2 ) = (831.9V/m) 2 [T(K)/300] 4 



(2) 



This electric field causes a temperature-dependent fre- 
quency shift of the atomic microwave clock transitions. 
It can be presented in the form (see, e.g. EI) 



Su/vq = p(T/To) 4 [1 + e(T/T ) 2 



(3) 



Here To is usually assumed to be room temperature (To 
300K). The frequency shift in a static electric field is 



5v = kE 2 . 



(4) 
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The coefficients k and (3 are reiated by 



= — (831.9V/m) 5 



(5) 



e is a smali correction due to frequency distribution QJ. 
In the present work we consider both terms in @ by 
calculating coefficients k, (3 and e. 

It is convenient to start from calculation of k by con- 
sidering an atom in static electric field. In the case when 
there is no other external electric field which sets pre- 
ferred direction the radiation shift can be expressed in 
terms of the scalar hyperfine polarizability of the atom. 
This corresponds to averaging over all possible directions 
of the electric field. The hyperfine polarizability is the 
difference of the atomic polarizabilities between different 
hyperfine structure states of the atom. The lowest-order 
effect is linear in the hyperfine interaction and quadratic 
in the electric field. Therefore, its value can be calculated 
using third-order perturbation theory (see, e.g. pr|) 



Se a = £ 



(a\V\n)(n\V\m)(m\V\a) 



(a\V\a)J2 



i\V\? 



(0) 



In our case the perturbation operator V is the sum of 
the hyperfine structure operator and the electric dipole 
operator 



V = H. 



hfs 



eE ■ r. 



The operator of the hyperfine interaction Hhfs is given 
by 



Hi 



hfs 



\e\ r x a 

— M • — 3 — , r> = max(r, r N ), (7) 



where a is the Dirac matrix, \i is the magnetic moment 
of the nucleus and rjy is nuclear radius. 

To get the effect of the electric field on the frequency 
of the hyperfine transition one needs to go through the 
following steps: 

• Substitute the perturbation operator V into equa- 
tion ©. 

• Keep only terms linear in Hhf s and quadratic in 
the electric field. 

• Apply equation JSJ to both components of the hy- 
perfine structure doublet. 

• Take the difference. 

The resulting expression for the frequency shift consists 
of three terms. The first two of them originate from 
the first term of equation |J?J|. In one of them the Hhfs 
operator is either on the left or right side of the expres- 
sion (these two terms are equal and can be combined 



together), and in the other the Hhfs operator is in the 
middle. The last term is due to change of the normaliza- 
tion of the wave function (second term of equation JBJ). 
It is proportional to the hyperfine structure of the ground 
state. 

After angular reduction these three terms become 



Sviias) = e 2 (E 2 ) 



21+1 
6 



A 



(ns||r||TOpj)(mpj||r||as) 



3 ) (^as 



(8) 
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(as\\r\\npj)A 



1/2) 



npj,mp] 



(mpj\\r\\as) 



n.m 



and 



^ ^ norm (^"^) — 

2/ P 2\2i + l, \ - \{as\\r\\m P] )\ 2 



(9) 



(10) 



Here 



c F = + + 1) - /(/ + 1) - ju + 1)] 

F' 

xj 1 ^ ] F' = \I-J\,I+J, 



A ns is the hfs constant of the ns state, A r 



is the off- 



diagonal hfs matrix element, I is nuclear spin, F = I + J, 
J is total electron momentum of the atom in the ground 
state ( J = 1/2 for atoms considered in present work), and 
j is total momentum of virtual p-states (j = 1/2,3/2). 
Summation goes over a complete set of ns, mpi/2 and 
mp 3 / 2 states. 

Expression JSJ does not include the s — d hfs matrix ele- 
ments while expression (|5J) does not include thepi/2~ P3/2 
hfs matrix elements. Test calculations show that the total 
contribution of the off-diagonal (in total momentum j) 
hfs matrix elements is of the order of 0.1% of the final 
answer and therefore can be neglected in present calcu- 
lations. 

Expressions © , © and i|10|) correspond to static limit 
when energy of thermal photon is zero. To take into 
account distribution one needs to make the following 
substitutions in terms JSJ) and ©: 



(TOPj||r||as) 1 
Ae^ ' 2 

and in term (|1C)|> 



(mpj||r||as) (mpj||r||as) 



Ae, 



Ae s 



At 2 



(Ae sp +uj) 2 (Aesp-cj) 2 



(11) 



(12) 



3 



where u> is the frequency of thermal photon. Integrating 
resulting expression with (yi and keeping only terms up 
to u! 2 (since ui <C Ae sp ) leads to expression of the form 
© in which first term is given by ©, © and 11U|) and 
parameter e in second term is given by 

Here index i replaces all indexes of summation in |J5j, 
© and 11UII . corresponds to terms in JSJ), k!f^ 

corresponds to 10, fc| 3 ' corresponds to (|1U|> and k = 

energy of the s — p transi- 
tion number i. If energies Ae sp j are in atomic units then 
A = 1.697 x 1CT 5 (atomic unit of energy is 27.211 eV = 
315773K). Lowest s — p transitions (e.g., 6s — 6p i/2 a nd 
6s — 6^3/2) strongly dominate in the summation (|13|) . 

III. CALCULATIONS 

In order to calculate frequency shift to the hfs transi- 
tions due to the electric field one needs to have a com- 
plete set of states and to have the energies, electric dipole 
transition amplitudes and hyperfine structure matrix el- 
ements corresponding to these states. It is possible to 
consider summation over the physical states and use ex- 
perimental data to perform the calculations. The low- 
est valence states for which experimental data is usually 
available dominate in the summation. Off-diagonal hfs 
matrix elements can be obtained to a high accuracy as 
the square root of the product of corresponding hfs con- 
stants: A m ^ n — y/A m A n (see, e.g. 0). However, the ac- 
curacy of this approach is limited by the need to include 
the tail contribution from highly excited states including 
states in the continuum. This contribution can be very 
significant and its calculation is not easier than the cal- 
culation of the whole sum. Also, for atoms like Yb + and 
Hg + available experimental data is insufficient to follow 
this path. 

Therefore, in the present work we use an ab initio ap- 
proach in which high accuracy is achieved by the inclu- 
sion of all important many-body and relativistic effects. 
We make only one exception toward the semiempirical 
approach. The frequency shift is dominated by the term 
(|10fl which is proportional to the hfs in the ground state. 
These hfs constants are known to very high accuracy from 
measurements for all atoms considered in the present 
work. It is natural to use experimental hfs constants 
in the dominating term to have more accurate results. 
Note however that the difference with complete ab initio 
calculations is small. It is also instructive to perform cal- 
culations of the hfs and atomic polarizabilities to demon- 
strate the accuracy of the method. 

Calculations start from the relativistic Hartree-Fock 
(RHF) method in the V* -1 approximation. This means 
that the initial RHF procedure is done for a closed-shell 



atomic core with the valence electron removed. After 
that, states of the external electron are calculated in the 
field of the frozen core. Correlations are included by 
means of the correlation potential method [2^] . We use 
two different approximations for the correlation poten- 
tial, E. First, we calculate it in the lowest, second-order 
of the many-body perturbation theory (MBPT) . We use 
notation E^ 2 ) for the corresponding correlation potential. 
Then we also include into E two classes of the higher- 
order terms: screening of the Coulomb interaction and 
hole-particle interaction (see, e.g. [29j for details). These 
two effects are included in all orders of the MBPT and 
the corresponding correlation potential is named 

To calculate E^ 2 ) we need a complete set of single- 
electron orbitals. We use the B-spline technique [30, |3ll 
to construct the basis. The orbitals are built as linear 
combinations of 50 B-splincs in a cavity of radius AOas- 
The coefficients are chosen from the condition that the 
orbitals are eigenstates of the RHF Hamiltonian Ho of the 
closed-shell core. The E^°°^ operator is calculated with 
the technique which combines solving equations for the 
Green functions (for the direct diagram) with the summa- 
tion over complete set of states (exchange diagram) (2t| . 

The correlation potential E is then used to build a 
new set of single-electron states, the so-called Brueckner 
orbitals. This set is to be used in the summation in equa- 
tions JHJl and (|10fl . Here again we use the B-splinc 
technique to build the basis. The procedure is very sim- 
ilar to the construction of the RHF B-spline basis. The 
only difference is that new orbitals are now the eigen- 
states of the Hq + E Hamiltonian, where E is either E^ 2 ^ 
or S(°°). 

We use the all-order correlation potential E^ 00 ) for Rb, 
Cs and Ba + . It has been demonstrated in a number of 
works (see, e.g. (2lJ IjU [3j| ) that inclusion of the screen- 
ing of Coulomb interaction and the hole-particle interac- 
tion leads to very accurate results for alkali atoms. How- 
ever, for other atoms with one external electron above 
closed shells these two higher-order effects are not domi- 
nating and their inclusion generally does not improve the 
results. Therefore, for the Yb + and Hg + ions we use only 
the second-order correlation potential E^ 2 ). 

Brueckner orbitals which correspond to the lowest va- 
lence states are good approximations to the real physi- 
cal states. Their quality can be checked by comparing 
experimental and theoretical energies. Moreover, their 
quality can be further improved by rescaling the corre- 
lation potential E to fit experimental energies exactly. 
We do this by replacing the Hq + E with the Ho + AE 
Hamiltonian in which the rescaling parameter A is chosen 
for each partial wave to fit the energy of the first valence 
state. The values of A are presented in Tabled Note that 
all values are very close to unity. This means that even 
without rescaling the accuracy is very good and only a 
small adjustment to the value of E is needed. Note also 
that since the rescaling procedure affects both energies 
and wave functions, it usually leads to improved values 
of the matrix elements of external fields. In fact, this is 
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TABLE I: Rescaling parameters for the £ operator which fit 
energies of the lowest s and p states of Rb, Cs, Ba, + , Yb + 
and Hg" 1 " 



TABLE II: Contributions to the hyperfine structure of the 
ground state of Rb, Cs, Ba, + , Yb + and Hg + (MHz); compar- 
ison with experiment. 



Atom 


a 


^1/2 


T)-i 


P3/2 


Atom 




Brueck 


RPA 


Str 


Norm 


Total 


Exp 


Rb 




0.868 


0.903 


0.906 


87 Rb 


5 s 


2888 


559 


-27 


-33 


3386 


3417° 


Rb 


Jj(oo) 


1.008 


0.974 


0.976 


133 Cs 


6* 


1957 


355 


-10 


-31 


2278 


2298 6 


Cs 


£(2) 


0.802 


0.848 


0.852 


137 Ba+ 


6* 


3509 


601 


-21 


-73 


4016 


4019 c 


Cs 


£(oo) 


0.985 


0.95 


0.95 


171 Yb+ 


6s 


11720 


1540 


-248 


-247 


12764 


12645(2) d 


Ba + 


±m 


0.782 


0.830 


0.833 


199 Hg + 


6s 


38490 


3873 


288 


-1483 


41169 


40507 6 


Ba+ 
Yb+ 


g(oo) 
£(3) 


0.988 
0.866 


0.963 
1.09 


0.967 
1.185 


"Reference 
''Reference 














Hg+ 


E< 2 ) 


0.805 


0.890 


0.926 


c References 
^Reference 


Tp and |3' 


71- 











e Reference 



a semiempirical method to include omitted higher-order 
correlation corrections. 

Matrix elements of the hfs and electric dipolc opera- 
tors are found by means of the time-dependent Hartree- 
Fock (TDHF) method [HHil- This method is equivalent 
to the well-known random-phase approximation (RPA). 
In the TDHF method, single-electron wave functions are 
presented in the form ij) — ?/>o + Sip, where ipo is unper- 
turbed wave function. It is an eigenstate of the RHF 
Hamiltonian Hq: (Hq — eo)V-'o = 0. Stp is the correc- 
tion due to external field. It can be found be solving the 
TDHF equation 



(14) 



where Se is the correction to the energy due to external 
field (Se = for the electric dipole operator), F is the op- 
erator of the external field (Hhf s or eE-r), and SV™^ 1 is 
the correction to the self-consistent potential of the core 
due to external field. The TDHF equations are solved 
self-consistently for all states in the core. Then matrix 
elements between any (core or valence) states n and m 
are given by 



(15) 



The best results are achieved when ip n and ip m are 
Brueckner orbitals calculated with rescaled correlation 
potential E. 

We use equation (|15|l for all hfs and electric dipole ma- 
trix elements in JSJ), 10, and (|10f> except for the ground 
state hfs matrix element in (|10f> where we use experimen- 
tal data. 

To check the accuracy of the calculations we perform 
calculations of the hfs in the ground state and of the 
static scalar polarizabilities of the atoms. Polarizabilities 
are given by the expression 



a (a) 



2^|(a|HK 



(16) 



which is very similar to the term (|10|l for the frequency 
shift. The most important difference is that the energy 



denominator is squared in term 1)10)) but not in 1)16)) • This 
means better convergence with respect to the summation 
over complete set of states for term 1)1011 than for 1)16)1 . 
Therefore, if good accuracy is achieved for polarizabil- 
ities, even better accuracy should be expected for the 
term HOI (see also Ref. 

However, the behavior of the other two terms, © and 
©, is very different and calculation of polarizabilities 
tells us little about accuracy for these terms. Therefore, 
we also perform detailed calculations of the hfs constants 
of the ground state. Inclusion of core polarization (second 
term in ()15)l ) involves summation over the complete set of 
states similar to what is needed for term (|SJ) . Comparing 
experimental and theoretical hfs is a good test of the 
accuracy of this term. 

In addition to term 1)15)) . we also include two smaller 
contributions to the hfs: structure radiation and the cor- 
rection due to the change of the normalization of the wave 
function. The structure radiation term can be presented 
in the form 



A str = (lpn\St\tpn) 



(17) 



where <5£ is the correction to the correlation potential 
due to external hfs field. The normalization term is 



<9S 

~A n (ip n \ — \ip n ) 



(18) 



where A n is given by 1)15)) with m = n. 

The results for hfs are presented in Table ITT1 Here col- 
umn marked as 'Brueck' corresponds to the (ip n \F\ip n ) 
matrix element. The column marked as RPA is the core 
polarization correction given by ("4> n \5V N ~ 1 \ r il) n ) , the no- 
tation 'Str' stands for structure radiation given by l)17)l . 
and 'Norm' is the renormalization contribution given by 
()18)l . In all cases ijjn is the Bruckner orbital correspond- 
ing to the ground state of the atom or ion, calculated 
with the rescaled correlation potential S. All-order E(°°) 
is used for Rb, Cs and Ba + . Second-order £( 2 ) is used for 
Yb + and Hg + . Comparing the final theoretical results 
with experiment shows that the theoretical accuracy is 
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TABLE III: Static polarizabilities a of Rb, Cs, Ba,+ , Yb 
and Hg + in different approximations (ajj). 



Atom 




E 




a c 6 


Total 


Other works 


S7 Rb 


5s 


E (2)c 


292.7 


9.1 


301.8 


329(23/ 








309.7 


9.1 


318.8 


293(6) s 






^(oo) e 


312.4 


9.1 


321.5 


318.6(6)'' 






AE (oo)d 


310.5 


9.2 


319.7 


318.5(6)' 


133 Cs 


6.s 


g(2)c 


343.8 


15.3 


359.1 


401.0(6) J 






A £(2)d 


383.5 


15.4 


399.0 


401. 5 ft 






E 


384.0 


15.5 


399.5 


400.4* 






AE (oo)d 


384.4 


15.5 


399.9 


400.6(1.0)' 


137 Ba+ 


6s 




104.1 


9.8 


113.8 








A £(2)d 


112.5 


9.9 


122.4 








^(oo) e 


112.8 


9.9 


122.7 








AE (oo)d 


112.7 


9.9 


122.7 




171 Yb+ 


6s 


g(2)c 


50.9 


6.2 


57.1 








AE (2)d 


55.4 


6.1 


61.5 




199 Hg+ 


6s 




10.5 


7.7 


18.2 








A £(2)d 


11.4 


7.6 


19.0 





"Polarizability due to valence electron. 
6 Polarizability of the core. 

C S( 2 ' is the second-order correlation potential. 

d R,cscaled E. See Table||]for the values of rescaling factors A. 

i s the all-order co rrelation potential. 
^Measurements, Ref. |2dj. 
9 Measurements, Ref. |2l|. 
^Calculations, Ref. u3- 
'Calculations, Ref. L'.l. 
^Measurements, Ref. [24j. 
^Calculations, Ref. |2fl. 
'Calculations, Ref. Illl . 

within 1% for all atoms except Hg + where it is 1.6%. If 
the structure radiation and normalization are neglected, 
accuracy for Rb and Cs remains within 1%, accuracy for 
Ba + becomes about 2% and accuracy for Yb + and Hg + 
becomes close to 5%. 

The results for polarizabilities, calculated in different 
approximations, are presented in Table HTT1 Pure ab ini- 
tio results obtained with E^ 00 ) and results obtained with 
rescaled correlation potential operators E^ 2 ) and E^ 00 ) 
are very close to each other and to other calculations 
and measurements. 



IV. RESULTS AND DISCUSSION 

Table HVI presents contributions of terms ©, 10 an d 
(|ll)fl into the total frequency shift of the hfs transitions 
for the ground states of 87 Rb, 133 Cs, 137 Ba+, 171 Yb+ 
and 199 Hg + calculated in different approximations. Term 
(f 1 Opt dominates in all cases, while term (JHJ) is small but 
still important at least for Rb, Cs and Ba + . Results ob- 
tained with E^ 2 ) and E^°°) differ significantly (up to 14% 
for Cs). However, after rescaling the results for both E^ 2 ) 



TABLE IV: Contribution of terms ©, ©, and JTUJ to the 

frequencies of the hyperfine transitions in the ground state of 
Rb, Cs, Ba+, Yb+ and Hg+ {Svq/E 2 x 1CT 10 Hz/(V/m) 2 ) in 
different approximations. 



A tAin 




V 
Zj 


1 S 1 

\ o 1 


\ d i 




Trvfal 
XOlcLl 


87pu 

Kb 


OS 


v( 2 ) a 


A C^C7 

-U.04o ( 


A A1 A 1 
U.U14 ( 


-u.ooyz 


-l.zUUo 






\V(2) b 


-U.ODOo 


U.U104 


-U.Doy4 


1 O/l AO 

-x.z4uy 






Zj 


A cc/in 
-U.004U 


A A1 


A 7AQ/1 
-U. /Uo4 


-l.ZDlo 






AZj ' 


-U.DOzU 


U.U104 


A SATO 


1 0/1Q7 

-l.z4o I 


133 nc 




v(2)a 


A G/1 1 O 

-U.y4iy 


A AO! A 

U.UzlU 


1 A799 


-l.yyol 






AZj 


-i.Uzoy 


A AOOQ 

u.uzzy 


-l.ZDoo 


9 9fiG7 






V^(oo) c 
Zj v ' 


-l.U14o 


A AOQO 

u.uzoz 


-l.Z i UD 


9 9^99 






\ * — > ( nri i h 

XT. K ' 


-1.0167 


0.0230 


-1.2695 


-2.2632 


137 Ba+ 


6s 


£(2)a 


-0.1027 


0.0034 


-0.1568 


-0.2561 






AE< 2 > 6 


-0.1095 


0.0036 


-0.1768 


-0.2827 








-0.1104 


0.0037 


-0.1778 


-0.2845 






xt, (oo)b 


-0.1104 


0.0037 


-0.1773 


-0.2841 


171 Yb+ 


6* 


£.(2)a 


-0.0672 


0.0009 


-0.0866 


-0.1529 








-0.0714 


0.0011 


-0.1003 


-0.1706 


199 Hg + 




£(2)a 


-0.0242 


0.0000 


-0.0296 


-0.0538 






xt (2) b 


-0.0263 


0.0000 


-0.0335 


-0.0598 



°S( 2 ) is the second-order correlation potential. 

''Rescaled E. See TableQJfor the values of rescaling factors A. 

c £(°°) is the all-order correlation potential. 



and E' 00 ) come within a fraction of a per cent of each 
other. Naturally, rescaling has a larger effect on results 
obtained with E^ 2 ). This means that the rescaling really 
imitates the effect of higher-order correlations and should 
lead to more accurate results. Comparing the results ob- 
tained with E^ 00 ), rescaled &°°> and rescaled E^ 2 ) gives 
a reasonable estimation of the accuracy of calculations. 
If we also combine this with the calculation of the hfs 
discussed above we can safely assume that the accuracy 
of the calculations for Rb, Cs and Ba + is on the level 
of 1%. Note that the frequency shift due to black body 
radiation can be a little larger (~ 1%) due to the effect 
of frequency distribution at finite temperature. 

For Yb + and Hg + we only have results with rescaled 
E< 2 ) and not rescaled t,(°°h They differ by about 11%. 
However, there are strong reasons to believe that the re- 
sults obtained with rescaling are more accurate. This is 
supported by calculations for Rb, Cs and Ba + as well as 
our experience with rescaling used in many earlier works. 
Therefore, the calculation of the hfs discussed above gives 
a more realistic estimation of the accuracy for Yb + and 
Hg + which is about 5%. 

Table llVl presents values of k (see formula |(2J). To 
obtain the frequency shift at finite temperature one needs 
to convert k into (3 using formula © and substitute (3 into 
equation ©. For accurate results one also needs to know 
the values of e. We calculate them using formula Ijl^jl in a 
very similar way as we calculate parameters k. Our final 
values of k, (3 and e are presented in Table [V] Parameter 
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TABLE V: Final results for the parameters k 
(KT 10 Hz/(V/m) 2 ), (3 (10~ 14 ) and e of the black-body 
radiation frequency shift for Rb, Cs, Ba, + , Yb + and Hg + . 



Atom 




k 


P 


c 


87 Rb 


5s 


-1.24(1) 


-1.26(1) 


0.011 


133 Cs 


6s 


-2.26(2) 


-1.70(2) 


0.013 


137 Ba+ 


6s 


-0.284(3) 


-0.245(2) 


0.004 


171 Yb+ 


6s 


-0.171(9) 


-0.094(5) 


0.002 


199 Hg+ 


6s 


-0.060(3) 


-0.0102(5) 


0.0005 



e for Cs was estimated in single-resonance approximation 
in Ref [lj and found to be 0.014. This value is in good 
agreement with our accurate calculations. 

The frequency shifts of some alkali metal atoms have 
been calculated and measured previously. We present 
previous results for the atoms and ions for which we per- 
form calculations in Table IVII together with our final re- 
sults. 

There is some disagreement for cesium. Our result is in 
good agreement with early measurements 0, 01 and ab 
initio calculations while recent measurements Id Ml 

and semiempirical calculations [a, llll Il2| give the value 
which is about 10% smaller. Less accurate measurements 
of Bauch and Schroder Q| cover both cases. We cannot 
comment on disagreement between experimental results. 
However, the source of disagreement between theoretical 
results seems to be clear. It comes from the contribution 
of the continuum spectrum to the summation over the 
complete set of states in term (JSJ. This term has off- 
diagonal hfs matrix elements between the ground state 
and excited states. Since the hfs interaction is localized 
over short distances (~ clq/Z) it emphasizes the contri- 
bution of states with high energies including states in the 
continuum (since ApAx > h, a small area of localization 
(Ax) allows high momentum (p) and thus high energy). 
In our calculations the contribution of states above 7p 
into term © is -0.35 x 10 -1 Hz/(V/m) 2 which is 15% of 
the total answer. 

In contrast, states above 7p contribute only about 
0.05% of the total value of term QlOfl. This is because 
the summation goes over the matrix elements of the elec- 
tric dipole operator which is large on large distances and 
thus suppresses the contribution of high-energy states. It 
is not surprising therefore that a semiempirical consider- 
ation, involving only discrete spectrum states, gives very 
good results for the atomic polarizabilities (see, e.g. [TT|). 
However, let us stress once more that the calculation of 
polarizabilities checks only term i|10|) and tells us very 
little about the accuracy of other two terms, © and JSJ. 

The contribution of the states above 7p is even more 
important for term @. Their contribution is about 30% 
of the total value of this term. However, the term itself 
is small and its accurate treatment is less important. 

In ab initio calculations by Lee et al [j| summation over 
complete set of states is reduced to solving of a radial 



TABLE VI: Electrostatic frequency shifts for the hyperfine 
transitions of Rb, Cs, Ba+, Yb+ and Hg+ (5u /E 2 x 10~ 10 
Hz/(V/m) 2 ) ; comparison with other calculations and mea- 
surements. 



Atom 


This 
work 


Other 

CdlC LLldLlOIlb 


Ref 


TV T nocn rnin ciri"(" t; 
IVItJeib HI clllcIlLo 


Ref 


87 Rb 


5s 


-1.24(1) 


-r.zuy4 


fal 

I9| 


-l.Zoyo ) 


I3| 


133 Cs 


6s 


9 OR(0\ 
-Z.ZUt Z ) 




rsi 


9 9Q^7^ 
-Z.zy^ ( J 


[91 








-2.2302 




-2.25(5) 


m 








-2.28 


[10] 


-2.17(26) 


f4] 








-1.97(9) 


M 


-2.271(4) 


f5] 








-2.06(1) 


[12] 


-1.89(12) 


f6] 








-2.268(8) 


[13] 


-2.03(4) 


\7] 


137 Ba+ 


(h 


-0.284(3) 


-0.27 


0] 






171 Yb+ 


6s 


-0.171(9) 










199 Hg+ 


6s 


-0.060(3) 


-0.058 









equation (equations of this type are often called Stern- 
heimer equation after one of the authors of this work). 
This approach does include the contribution of the con- 
tinuum spectrum and the result is in very good agreement 
with ours (see Table fvT)l . 

In other ab initio calculations by Pal'chikov et al [To| 
summation is done via Green functions. This corresponds 
to summation over the complete set of states and does 
include the continuum spectrum. Again, the result is in 
very good agreement with other ab initio calculations (0 
and the present work). 

Recent calculations by Beloy et al 01 applied a mixed 
approach, with extensive use of experimental data for 
lower cesium states and ab initio summation over higher 
states including continuum. The result is in good agree- 
ment with fully ab initio calculations. 

In contrast, analysis performed in H,El^2 i s limited 
to discrete spectrum. Adding —0.34 x 10 _1 Hz/(V/m) 2 
(which is total tail contribution from all three terms JSJ, 
© and HI Up found in our calculation) to the results of 
Feitchner et al and Micalizio et al brings them 
to excellent agreement with ab initio calculations. The 
same modification of the result by Ulzega et al 12] makes 
it a little bit too large but still closer to other results than 
without the tail contribution. 



V. CONCLUSION 

We have performed calculations of the frequency shift 
of the ground state hyperfine transition for several atoms 
and ions caused by a static electric field which can be 
used to evaluate the effect of blackbody radiation on the 
frequency of the microwave atomic clock transitions. The 
size of this shift is comparable to the current error in the 
measurements of the energy shift caused by variation of 
a and so needs to be taken into account in laboratory 
measurements placing limits upon a variation. 
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Detailed analysis of the calculations for cesium reveal 
the source of disagreement between different theoretical 
approaches. This seems to be contribution of the con- 
tinuum spectrum into summation over complete set of 
states which was neglected in semiempirical calculations. 
Restoring the tail contribution in works where it was ne- 
glected brings all theoretical results in good agreement 
with each other. 
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